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Abstract: Motivated by recent developments in superstring theory in the cosmological 
context, we examine a field theory which contains string networks with 3-way junctions. 
We perform numerical simulations of this model, identify the length scales of the network 
that forms, and provide evidence that the length scales tend towards a scaling regime, 
growing in proportion to time. We infer that the presence of junctions does not in itself 
cause a superstring network to dominate the energy density of the early Universe. 



Keywords: 



Superstrings, cosmic strings, cosmology . 



1. Introduction 



Type IIB string theory, after compactification to 3+1 dimensions, has a spectrum of 1- 
dimensional objects in the form of solitonic D-strings and fundamental or F-strings which 
are of great current interest in cosmology Q H|. Such objects offer a much richer form 
of dynamics than the more commonly studied Nielson-Olesen string of the Abelian Higgs 
model. Individually, the F and D strings are ^-BPS objects, but they each break a different 
half of the super symmetry. When combined in a bound state one may think that this would 
further break the super symmetry, but in fact the bound states are still 5-BPS and give the 
familiar BPS tension for bound states, 



VFVp 2 + q 2 /gl (i-i) 



Where we have p F-strings and q D-strings, \ip is the effective fundamental string tension 
after compactification and g s is the string coupling. Given that stable bound states of 
vortices are present we should expect that junctions will also be in the spectrum, where 
two different types of string come together at a point and form a bound state leading away 
from that point. Such junctions not only exist, but are j-BPS states ]|, |J. In fact, it is 
possible to construct whole webs of these junctions and still preserve supersymmetry ||. 
This is not the whole story as far as stability is concerned, for example we know that strings 
can break by ending on a D-brane. In Q the authors give a comprehensive account of the 
mechanisms which can lead to cosmic superstring breaking, concluding that it is possible 
for stable cosmic-superstrings to exist. Knowing that such strings can exist is not the same 
as saying they do: a mechanism is required to generate them. An attractive scenario is 
provided by the annihilation of a D3-D3 pair terminating a period of brane inflation, which 
results in the formation of a networks of (p, q)-strings || [jj. 

In warped compactifications /j,f can take any value below the 10D string tension Tp, 
but in brane-antibrane inflation models the tension is related to the inflation scale, with 
G/i ~ 1CP 6 . String networks surviving until today promise observable signals such as 
Cosmic Microwave Background (CMB) fluctuations and gravitational lensing. Indeed, 
much excitement was generated by the observation of an unusual double galaxy, thought 



to be a candidate for the first observation of a string lens ||, ||, 10 1, but unfortunately ruled 
out by subsequent high-resolution observations fpTjJl . 

The dynamics of networks even of ordinary cosmic strings in are not well understood. 



In the conventional scenario [12, |l3j energy is lost from the network into loops, which oscil- 
late and slowly decay into gravitational radiation. However, when an underlying (classical) 
field theory with string solutions such as the Abelian Higgs model is simulated the energy 
goes into gauge and Higgs radiation via some kind of non-perturbative process acting at the 
string width scale [14, |Tq ]. Whatever the energy loss mechanism, there is broad agreement 



that a network of cosmic strings formed in the early Universe will settle into a self-similar 
or so-called scaling regime, in which the gross length scales such as the average inter-string 
distance and the average curvature radius increase in proportion to time. 

Scaling is an extremely important property, for it ensures that the strings contribute a 
constant fraction of the energy density and therefore do not come to dominate the Universe 
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(or disappear to a negligible density). The string density p s is approximately given by the 
string tension multiplied by the length per unit volume, neglecting kinetic energy. On 
dimensional grounds, we can define a length scale £ such that p s = /u/£ 2 , and scaling 
means that £~torp~^/i 2 . Hence dividing by the critical density p c oc 1/Gt 2 of a 
Friedmann universe, we see that Q s ~ G\i. 

However, once junctions are added to the network, it is not clear that scaling will be 
maintained, which would be cosmologically disastrous for the new cosmic string scenarios. 
A "frozen" string network would act like a fluid with equation of state p = — p/3 [16 



17] which cannot be reconciled with current cosmological data. It is therefore extremely 
important to investigate the dynamics of such networks. 

The first simulations of string networks with junctions Jl^, |l^j used very simple models 
through which it was hoped the dynamics of the field theory could be captured, and gave 
respectively negative and positive indications of scaling. More recent attempts at modelling 
of string networks with junctions |2(| give indications that scaling is allowed. Where there is 
an underlying field theory, it is to be expected that direct simulations of the field equations 
will give more reliable results. Here again results are mixed. There is evidence from 
numerical simulations of domain wall networks in 2+1 dimensions both of departures from 



naive (linear) scaling [21, 22] and consistency with it The issue has been addressed 



recently in 3+1 dimensions p4| ], reviving earlier work of pl|] , where the non- linear sigma 
model (NLerM) approximation of a field theory containing multiple vortices and junctions 
was shown to scale with the background fluid. The NLuM approximation greatly reduces 
the number of dynamical variables and makes it possible to perform large-scale simulations. 
However, the model is somewhat peculiar in that the c-model space is not in fact a manifold, 
and in any case misses the microphysics associated with the strings. 

Here we perform for the first time full field theory simulations for a model of a string 
network with junctions, specifically Z3 vortices, and show that the evolution is consistent 
with a late-time scaling regime. While our model does not fully capture the dynamics of 
(p, g)-strings, we can conclude that the presence of junctions is not in itself inconsistent 
with scaling and, by extension, that (p, (^-string networks are viable for sufficiently low 

PF- 

2. The model 

A field theory with spontaneous symmetry breaking from Lagrangian symmetry group G 
to vacuum symmetry group H has string solutions if the manifold A4 ~ G/H (which may 
not coincide with the manifold of minima of the tree- level potential [26]) is not simply 



connected, or its first homotopy group is non-trivial: tti(A4) 7^ 0. If G is simply connected 
this can be ensured by arranging for H to be disconnected: ttq(H) 7^ 0. 

Disregarding the Abelian symmetry-breaking G=U(1), H=l, the smallest simply- 
connected compact Lie group is SU(2), which suggests its use as the symmetry of the 
Lagrangian of a field theory. The most trivial discrete group is the cyclic group Z n and so 
we are led to the conclusion that the simplest model with string junctions in a non- Abelian 
field theory would have a vacuum manifold SU(2)/Zat. To obtain a vacuum manifold 
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SU(2)/Zat we use an adjoint scalar, <J> = <& a T a and a complex, (2j + l)-component vector 



(ft, where j is the spin of the representation, which transform as 

$ U$U~ l (2.1) 

<f>-^U<f> (2.2) 

The action we use is given by 

S = j d 4 x(T- V) (2.3) 

where the kinetic energy density takes the standard form 

T = -d^d^cft - Tr^*^*) (2.4) 
and the potential is broken down into three pieces V = Vi + V2 + V3, with 

Vi = Ai(Tr(cI> 2 ) - rj 2 ) 2 + \\ 2 {(ftH - v 2 ) 2 , (2.5) 

V 2 = -A 3 # f $0 - A 4 (^ T C T + - A 5 (0 T C T $0 + ^t$c*0*), (2.6) 

V 3 = \ 6 \4> T C T (ft\ 2 + A 7 0VTr(cD 2 ), (2.7) 



where C is the charge conjugation matrix (see Appendix |A|). The first potential, V\ is 
used such that the fields aquire a non-zero vacuum expectation value (vev). The second 
potential is present in order to get the two fields to interact, and explicitly break the residual 
U(l) symmetry of V\. The third part is included for completeness, and is not necessary 
to obtain the desired symmetry-breaking. We accordingly set \q and A7 to zero in the 
following. Note that if j is half-integer, the charge conjugation matrix is skew-symmetric 
and so the A4 and A6 terms vanish. 

3. Vortex solutions 

Our (2j + l)-dimensional SU(2) basis T a is chosen so that 

Tr(T a T b ) = V, T 3 = ^-diag(j, j - 1, . . . , -j + 1, -j), (3.1) 

with N 2 = § (2 j + + 1) (see Appendix [A]). 

We can understand the presence of vortices in the model by considering the vacuum 
structure. When the fields are in the vacuum, we can use two of the SU(2) rotations to get 
$ = hr/T 3 , where h is a dimensionless constant. Taking (ft = v (a(j), a (j-i)> •••> a (-j+i)> a {-j)) 
we find that the first term in V2, 

-A 3 # t $0 = -^A 3 ^ 2 (j>(i)l 2 + U ~ ^l^i-i)! 2 + •••) ( 3 - 2 ) 

is minimized for a given (ft' (ft when aj is the only non- vanishing component. We can use the 
final SU(2) rotation to ensure that aj is real giving the following ansatz for the vacuum 

$ = hr]T 3 (3.3) 
(ftl = v(a (j) ,0,0,...) (3.4) 
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Figure 1: A representation of the vortex profile functions \(j)\ 2 and Tr<I> 2 , with parameters given 



in Eq. (4.4) 



Having established the form of the vacuum, we are in a position to find its residual sym- 
metry. Performing a rotation with T 3 = NjT 3 , Uq = exp(i#T 3 ) we find 



$0 

exp(ij6)4> 



(3.5) 
(3.6) 



so we see that the vacuum is unchanged for 9 = 2irn/j, which gives the discrete structure 
Z^r we were aiming for, with N = 2j. As we have a Z^r structure for the vortices then we 
see that there must be junctions. To understand this we consider the case for Z3, which is 
the model we shall actually be simulating. In two space dimensions, a finite energy field 
configurations representing three vortices together can be continuously deformed to the 
vortex-free vacuum, while maintaining finite energy. Similarly, two anti-vortices can be 
continuously deformed to a vortex. If we make the deformation parameter to depend on 
a third space dimension, and we realised that there will be junctions where the transition 
takes place, presumably localised in order to reduce the total energy. One can also view 
the vortices as the bound-state of the anti-vortices of the other two, as for (p, q) strings. 

Another way to understand junctions in our model is by noting that we have an adjoint 
SU(2) scalar field, and therefore when v 2 < in the potential (|2.5| ) the unbroken symmetry 
expands to U(l), and there will be monopole solutions [^7|]. An expectation value for the 
field (j) breaks this U(l) down to Z^v, and confines the monopole flux to N strings. Hence 
the monopoles become points where three vortices come together to form a junction. Note 
that in the gauged Z2 case the monopoles are "beads" on strings [28|. 

This is far from being the only way to obtain vortices with junctions from a field theory. 
Indeed, the Abelian symmetry breaking U(l)— >1 has a spectrum of vortices labelled by the 
integers, and stable bound states of vortices exist both in the pure scalar theory, and in 
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the gauge theory when the gauge boson is lighter than the Higgs boson (corresponding to 
a Type II superconductor). Recent work has shown that string networks with junctions 
can form in this theory (2t|, but not whether the networks scale. 

A model closer to the (p, g)-string, which has two separate U(l) symmetries corre- 
sponding to the electric flux of the F-string and the magnetic flux of the D-string, has 



recently be proposed by one of the authors [30]. There are two U(l) symmetries in this 



model, but both kinds of string have magnetic fluxes, and the characteristic string tension 
formula (|l]l]) is not exactly obeyed. 

A more complex model based on a gauged U(N) with a global SU(N) symmetry has 



also been put forward [31]. The global SU(N) symmetry can be lifted in such a way as to 
give string junctions, and the strings have the very interesting feature that it is possible 
for different kinds to pass through one another without reconnecting. 

4. Evolution 

Ideally we would like to perform the simulations in the background of some standard 
cosmology, typically a radiation or matter dominated Friedmann model. In practice it 
turns out that for field theories with an explicit mass scale this is problematic. In order to 
properly simulate the vortices one must choose a sufficiently small lattice spacing to resolve 
their core throughout the whole simulation. In an expanding background, the physical 
separation between lattice sites increases with time, so the vortices effectively "fall through 
the lattice" after some period of time. The other effect of expansion is Hubble damping, 
which is easy to simulate. In our simulations we wish to mimic the effect of expansion in 
a radiation dominated Universe by first transforming to conformal time, r, the equations 
of motion for a real scalar with potential V then become 

ip " + 2 ^ ip '-V 2 ^ = -a 2a ^-V, (4.1) 
a oip 

where the prime denotes differentiation with respect to r, a(r) is the scale factor, V 2 is the 
co-moving Laplacian, and a is a constant which should be 1. However, as explained above 
this causes problems at late time with the vortices not being resolved on the lattice. Instead 
we adopt a tactic suggested in [^] , which is to reduce the rate at which the strings shrink 
by reducing the size of the parameter a. In our simulations we take a = 0, which results 
in vortices of constant comoving size, we also evolve with the Hubble damping appropriate 



to a radiation dominated Universe. According to 32], this approximation has little effect 
on the dynamics of the vortices themselves, at least in the U(l) model they considered. 

String networks are thought to quickly evolve towards a statistically simple self-similar 
scaling solution, where all the relevant length scales in the problem evolve proportional 
to time. What we aim to determine here is whether non-Abelian networks of the type 
constructed above follow the scaling rule. The most elementary scale of the problem is 
the average separation of strings, £, and this can be constructed by knowing the average 
comoving length L of string per unit co-moving volume V as 

l/e s = L/V. (4.2) 
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Scaling would mean that £ s oc r. 

Determining the length of string in this model in a lattice simulation is not as simple 
as for the Abelian model with a simple complex scalar field, where the field vanishes in the 
core of the string in the continuum. In that case we can identify plaquettes through which 
a string passes by use of the geodesic rule, whereby the phase between two neighbouring 
lattice points is assumed to interpolate between the two values using the shortest path. No 
such convenient local indicator of the string position is guaranteed to exist on a lattice, so 
we instead use the profile functions of the vortices, Fig. |l], noting that \(p\ is at a minimum 
(but not zero) at the centre. We therefore identify the length of a string by the volume of 
the region where \4>\ 2 < |</> c | 2 , divided by the area of the static solution with |</>| 2 < \4> c \ 2 , 
for some cutoff \4> c \ 2 - In the early stages of the simulation we would not expect this to 
be a good measure due to the large oscillations of the field, however, as the simulations 
proceeds these oscillations decay away thanks to Hubble damping, leaving a cleaner signal. 
In practice we took three different cutoffs, finding that there was no significant difference 
between them, lending credence to this as a reliable measure of string length. 

For the junctions we can construct a similar length scale using dimensional analysis, 

Vd = N/V, (4.3) 

and this corresponds to the typical distance between junctions. The subscript m is used 
because the junction charge corresponds to the monopole charge of the adjoint scalar 
N/V is then the number density of monopoles. Unlike the string charge, we can evaluate 
the integer monopole charge on a lattice, and in the simulations we use this to measure N. 

We now come to the details of the simulations. The equations were solved using a 
simple 3-point discretisation of the Laplacian and a leapfrog time evolution, with the code 
parallelised using the LATfield library for lattice field evolution We used a 544 3 box 
with the scalar potential parameters 

Ax = 3/4, A 2 = l, X 3 = N j /j, n = u = \ (4.4) 

and the rest either irrelevant or vanishing, as discussed above. The lattice spacing was 
taken as h x = 0.4 and the time step given by ht = Q.3h x . The simulations were run up 
until the time it takes for light to cross half of the box, using a radiation background as 
discussed earlier. For initial conditions we used the field configuration which resulted from 
setting the fields to their vacuum value with random orientations, then damped to remove 
the large amount of radiation. The damping lasted for 200 time steps and the damping 
factor during this period was V = 2.0, where the damping is included as <f>" + T(f>' + the 
half-box crossing time includes this damping period. 

We performed nine realizations with a typical snapshot as shown in Fig. where we can 
clearly make out the prescence of vortices and junctions. The result of these simulations 
is given in Fig. ||, showing the mean evolution of the length scales £ s and £ m , with 1-a 
statistical errors either side. From this plot we see the linear behaviour of both length 
scales indicating scaling behaviour. To evaluate the gradients, £, we present Fig. ||| showing 
that at late times the gradient is consistent with being a constant, implying scaling for £ s 
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Figure 2: a snapshot of a run, showing the three-point vertices in a run using the parameters of 



and £ m , with scaling value extracted from the last half of the simulation 



Cs = (0.4±0.1)r, 
■ m = (0.5 ± 0.2)r. 



(4.5) 
(4.6) 



5. Conclusion 

We have carried out numerical simulations of a classical field theory containing strings 
with junctions, in order to test whether the presence of junctions in the network affects 
scaling. Our model is based on the global symmetry breaking SU(2)— >U(1)— ► Zjv, and 
the junctions can be thought of as global monopoles connected by global strings. We find 
good evidence, displayed in Fig. |3|, that the string network length scale £ s (defined in Eq. 
|4.2j ) scales in proportion to time with coefficient 0.4 ± 0.1 and the junction length scale £ m 
(defined in fl4.3| )) scales with coefficient 0.5 ± 0.2. 

While our model was motivated by an attempt to model (p, q) string networks, it fails 
to capture some important features:the spectrum of string tensions differs, with the field 
theory giving equal tension to all vortices; the reconnection probability is close to or equal 
to unity. Our conclusions are that it is perfectly possible for networks with junctions to 
scale, backing up the results for the NLaM approximation in [EJ], and that they do not 
in themselves constitute a cosmological disaster. However, we must wait for better models 
before drawing quantitative conclusions for (p, q) string networks. 
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time 

Figure 3: Top: network length scales £ s (solid) and £ m (dashed), against conformal time r, 
averaged over nine realisations. Bottom: the time derivatives of the network length scales, again 
averaged over nine realisations. 
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A. SU(2) conventions 

To construct the spin-j representation of SU(2) we begin with a basis for the algebra J±, 
J3 normalised in the usual way 

[J 3 ,J±] = ±J±, (1.1) 
[J+,J_]=2J 3 . (1.2) 

We identify the Casimir J 2 = J_ J+ + J3(J-s + 1) and construct normalized eigenvectors of 
J 2 and J3, 

J 3 \P,m) = m\P,m) (1.3) 
J 2 \P,m) = p\p,m) (1.4) 
{P,m\P,m) = 1 (1.5) 
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By considering the highest and lowest weights one finds the following matrix components 
for the generators T 1 



^m,m = (j,m\J + \j,m) = - rh){j + rh + l)5 m ,f. 



- m,m 
r3 



U,m\J-\j,rh) = y/{j + fh){j -fh + l)5 m ,f. 
(j,m\J 3 \j,rh) = rhSm^rh 



(1.6) 
(1.7) 
(1.8) 



with (3 2 = j(j + 1) and -j < m,rh < j. Defining T 1 = (f+ + f~)/2, f 2 = (f+ - f~)/2i 
we have 



a 

In our model we rescale by Nj, where 



(1.9) 



2 
3 ( 



]V/ = ~(2.7 + l)i(j + l) 



to get 



rna _ rpa 

Ti(T a T h ) = ^S ab . 



(1.10) 

(1.11) 

(1.12) 



In the model we also use the charge conjugation matrix C which is the intertwiner between 
the SU(2) representations T a and — T a *, 



So, by defining SU(2) group elements as U = exp(i6 a T a ) we have 

CU = U*C. 



(1.13) 



(1.14) 



In our basis we take 



C 



( ... l\ 
..0-10 
..100 



V 



(1.15) 
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